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We show that the rate of separation of two phases of dif- 
ferent densities (e.g. gas and solid) can be radically altered 
by the presence of a metastable intermediate phase (e.g. liq- 
uid). Within a Cahn-Hilliard theory we study the growth in 
one dimension of a solid droplet from a supersaturated gas. 
A moving interface between solid and gas phases (say) can, 
for sufficient (transient) supersaturation, unbind into two in- 
terfaces separated by a slab of metastable liquid phase. We 
investigate the criteria for unbinding, and show that it may 
strongly impede the growth of the solid phase. 
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The influence of metastable states on the dynamics of 
phase transformations is an important issue in materials 
physics. Their role has long been acknowledged in met- 
allurgy, where the rate of phase transformation is often 
limited by conduction of (latent) heat |l|,||. Indeed, Ost- 
wald's empirical 'rule of stages' || asserts that the trans- 
formation from one stable phase to another proceeds via 
all metastable intermediates in turn. In complex fluid 
systems such as colloidal suspensions, particle diffusion, 
not heat diffusion, is often rate-limiting Here too, 
a strong influence of metastable phases on ordering dy- 
namics has been suggested, e.g. in the phase ordering 
kinetics of polymer liquid crystals ^ , random-coil poly- 
mers H , proteins M and colloid-polymer mixtures . 
In this Letter, we propose a simple model to account for 
the role of metastable states in phase ordering limited by 
particle diffusion (described by a conserved order param- 
eter; for the nonconserved case, e.g. nematic order, see 
Refs. ||). 

To establish ideas, consider first a homogeneous fluid 
of a simple substance which is quenched to a temper- 
ature below its triple point. It will separate into two 
coexisting phases, the gas and the solid, whose densi- 
ties are given by the construction in Fig. However, in 
this temperature range, the low density (fluid) branch of 
the free energy has an additional minimum, representing 
the metastable liquid phase, whose presence can interfere 
with the phase separation process. Free energy curves of 
precisely this form can be realized in several complex fluid 
systems, notably mixtures of spherical colloids and much 
smaller polymers ]TT| . Experiments Q show that, when 
such a mixture is prepared with a composition at which 
a metastable minimum is present JTT| it does not sepa- 
rate quickly into dilute colloidal fluid and dense colloidal 
crystal (the predicted equilibrium phases). Instead, af- 



ter an initial latency period (consistent with nucleation) 
dense but non-crystalline droplets separate out, forming 
an amorphous sediment which begins to crystallize only 
on a much longer time-scale Jl2| . 

We can explain these qualitatively robust results 
within a simple, one-dimensional continuum model if we 
take due account of the metastable 'liquid' minimum in 
the free energy curve. (We now use 'gas, liquid, solid' 
as generic labels for three phases of increasing densi- 
ties (|l2|.) Details of our results will be given elsewhere 



Consider a system with a conserved order parameter 
(a concentration or density (/> Q) at equilibrium. At co- 
existence, there is an interface between two phases ('gas' 
and 'solid') whose densities (<j) g and (f> s ) are found from 
the bulk free energy density /(</>), as in Fig. [|, by the 
double-tangent construction which ensures equality 
of the chemical potential \x and the (osmotic) pressure P 
in the two phases, (/i is the slope of the tangent and P 
its negative intercept.) In minimizing the free energy, we 
seek the lowest possible double tangent, and may safely 
ignore local minima such as that labeled 'liquid' in Fig. |]. 
Across the interface, <fi changes smoothly (as in Fig. ||a), 
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FIG. 1. Bulk free energy density / versus concentration 
(/>, with double tangents for equilibrium binodals (solid line), 
and metastable binodals (dashed lines). The supersaturation 
a for an initial homogeneous state (•) is denned as shown. 

passing through values for which the \x and P calculated 
from f(4>) are "wrong". However, the curvature of the 
concentration, V 2 </>, contributes to /i and P so that both 
are constant for the smooth interface which exists in the 
thermodynamic ground state. 

Similar interfaces form in the early stages of phase sep- 
aration, between regions of different cj) jlB ], In the so- 
called "growth stages" <p(x_,t) evolves primarily by mo- 
tion of these interfaces |16| . In studying this motion, 



1 



metastable minima are typically ignored but, as we now 
show, they can play a crucial role. 

For simplicity we consider only one spatial dimension 
(perpendicular to an interface) and study the evolution 
of the interfacial density profile. The Cahn-Hilliard equa- 
tion llTtl for this time evolution is 



Conc.(|) 



dt dx \ dx \ d(f> dx 2 J J 



(1) 



where the mobility T and curvature constant K are phe- 
nomenological parameters; we are interested in f{4>) of 
the form shown in Fig. [y. Eq. ([j]) is derived by equat- 
ing dcj)/dt to the divergence of a current J = -rV/j 
induced by a chemical potential gradient. The chemical 
potential is defined in terms of a functional derivative, 
|U = 6F[4>]/6<j>, where F[cj>] is the integral of the local free 
energy density, f{4>) + K(\7<j)) 2 /2. 

Consider first an ordinary, smooth solid-gas interface 
as shown in Fig. ||a. If the solid region is to grow, the in- 
terface must move to the right, and this requires a chem- 
ical potential gradient in the supersaturated gas phase. 
Steady-state solutions of Eq. ([j]) (a uniformly moving 
interface) are thus unphysical in far field, since a non- 
zero gradient would have to extend to infinity. However, 
such solutions do shed light on the local dynamics of the 
interface since, during the growth stage p8[ , local re- 
arrangements are much quicker than changes in the far 
diffusion field driving the interfacial motion. Hence the 
interface experiences a quasi-constant flux from the su- 
persaturated region; its local dynamics are quasi-steady- 
statc Q . Steady-state solutions of Eq. (|l]) are calculable 
to arbitrary accuracy for any form of f(<j>) and r(c/f)), and 
exactly for piecewise-quadratic f(<p) and piecewise con- 
stant r Jl3|]. (Constant K is assumed.) In the latter case, 
linear stability of an ordinary moving interface (Fig. ||a) 
can be established rigorously Despite this, we now 
explore an alternative, 'split' mode of interfacial motion. 

Since the double-tangent construction identifies states 
of equal chemical potential and pressure, an interface 
may form between points on the metastable binodal, such 
as (<Pa,<I>b) in Fig. |l|. So, a split profile such as that 
shown in Fig. ||b, with well-separated gas-liquid (g/1) 
and liquid-solid (1/s) parts, has only gentle gradients of 
\i, and is therefore long-lived. Physically, a moving gas- 
solid (g/s) interface (Fig. ||a) is locally stable against this 
splitting because, if the flux of material from the 
right onto the g/s interface is increased (so as to advance 
the base more than the top, splitting the interface) the 
curvature, V 2 (/>, is perturbed, setting up currents tend- 
ing to restore the profile. The g/1 and 1/s interfaces can 
be described as curvature-bound. This negative feedback 
mechanism does not operate for the well-split interface 
of Fig. |^b. In the metastable liquid region, curvature is 
small, and the Cahn-Hilliard equation reduces to the dif- 
fusion equation Jl5|. Material flows to the 1/s interface 
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FIG. 2. (a) A solid-gas interface. Note the inflection near 
the density of the 'liquid minimum' in f(4>). (b) A split in- 
terface: the inflection has become a slab. Solid, metastable 
liquid and gas regions are separated by 'curvature-unbound' 
interfaces at x = xi(t) and X2(t), between binodal densities 
((I>a,(I>b) and (4>c,4>d) respectively (see Fig. 0). Fluxes ja, 
Jb, jc flow in and out of the interfaces. In the gas phase, the 
supersaturation relative to the metastable binodal 4>a, is a'. 

from the g/1 interface down the intervening concentration 
gradient, whose sign follows directly from the metastabil- 
ity of the liquid, as can be seen from the double tangent 
constructions in Fig. If the flux of material onto the 
g/1 interface is sufficient to increase its separation from 
the solid droplet, the concentration gradient in the mid- 
dle region, and hence the restoring flux, will be reduced. 
So a positive feedback mechanism operates; the split in- 
terface is then curvature-unbound. 

We have observed curvature-unbinding in numerical 
solutions of Eq. ([!]), using f(4>) of the form sketched 
in Fig. 0. A time-sequence is shown in Fig. |^; the ini- 
tial condition is a sharp step between a solid at density 
4> s and a supersaturated state (</> = 4> g + a) within the 
gas (rather than liquid) well of the free energy curve. 
Below some critical supersaturation (discussed below), 
such curvature-unbinding does not occur, and the inter- 
face propagates in the conventional (unsplit) mode. 

We now analyse the dynamics of a pair of curvature- 
unbound interfaces. On scales much larger than the char- 
acteristic interfacial width y/K/f", the g/1 and 1/s in- 
terfaces can be represented by sharp steps at positions 
xi(t) and X2(t) respectively (Fig. ^d). Concentrations 
at the interfaces will be fixed at the metastable binodal 
values (Fig. |l), so long as interfacial velocities are small 
(x\,x 2 <C \ f" 3 / 2 K~ 1 / 2 T). In the gas and liquid regions, 
evolution is purely diffusive, with diffusion constants D g 
and Di respectively; D is set to zero in the solid phase. 
Conservation of material requires (4>b — 4>a)x\ = ja — 3b 
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FIG. 3. Superimposed snap-shots at equal time intervals, 
from numerical solution of Eq. (^), using f(4>) as sketched in 
Fig. |lj Initial condition is a sharp step between <f> s and 4> g + a. 
A region of metastable liquid forms as the two parts of the 
interface "curvature- unbind" . Inset: Wide- frame view of the 
last time-slice. 
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FIG. 4. Size of the metastable region x\ — X2 against time 
t. Solutions of finite lifetime (dotted lines) for a' /a c = 0.58, 
0.70 and 0.75. The marginal solution a = a c (dashed line). 
Diffusively unbound solutions (solid lines) for a' /a c — 1.1 and 
3.2. In each case x 1 — X2 = 1 at t = 0. 



and {4>d — 4>c)i2 = jc, where the fluxes ja, 3 b an d jc 
are defined in Fig. ||b. 

The control parameter is now the supersaturation rel- 
ative to the metastable binodal a' = a — ((f) a — <ftg) ( see 
Fig. |IJ). Two further approximations are made, whose 
adequacy has been checked by numerical solution of the 
full equations fM| , The first is that the diffusive field in 
the gas phase is unperturbed by the motion of the steps, 
leading (lj] to the equation ja ~ o' sj D g /nt. The second 
is that the gradient in the liquid region is approximately 
uniform, leading to j B « jc ~ Di (0b - <f>c)/{xi - x 2 ). 
These yield tractable equations of motion for the size of 
the metastable liquid region: 
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and for the size of the solid droplet: 
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Here 



<p A ) 2 /(D g o-' ) and a c = 



D t /D g } 1/2 . 
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For < a 1 < cr c , solutions of Eq. (||) collapse to zero 
size in a finite time, as might be expected of a metastable 
phase. We describe these interfaces as diffusively bound. 
(This is a much "looser" binding than the curvature bind- 
ing discussed above.) After collapse of the metastable 
slab, the interfaces will curvature-recombine, and the 
solid will resume the ordinary (unsplit) mode of growth. 
In contrast, for a' > a c solutions of Eq. (ph have open 



trajectories. The interfaces are diffusively unbound, and 
the metastable liquid region grows, in principle, for ever. 
Typical trajectories are shown in Fig. ||. 

In the unbound case, to reach the solid, material must 
diffuse through the liquid phase (down a gradient which 
becomes flat at the triple point o~ c — > 0); the growth 
rate of the solid region is thereby suppressed, by a factor 
which we find to be of order (a c /a') 2 , compared to the 
case of a curvature-bound interface [Q. In this way, 
solid growth from a supersaturated gas can be strongly 
arrested by a weakly metastable liquid phase. Fig. || 
shows the size of the solid region x 2 at a given late time 
to, as a function of supersaturation a. For a' > o~ c , the 
solution of Eq. (||) is shown. For a' < a c , the metastable 
region collapses at some early time, after which ja ~ 
o-^/Dg/nt. 

This scenario of diffusive unbinding of interfaces and 
suppression of crystal growth relies on the initial forma- 
tion of curvature-unbound interfaces during the nucle- 
ation stage. Whereas the critical supersaturation a c for 
diffusive unbinding is independent of initial conditions, 
that for curvature-unbinding defines a region, in the 
space of initial configurations, with a nontrivial bound- 
ary. In our (non-exhaustive) numerical study of evolution 
of a step function, curvature-unbinding appeared to oc- 
cur whenever a' > cr c ; if so, for these initial states at 
least, the unbinding is diffusively controlled. 

Real systems are not infinite and hence the ambient 
supersaturation must eventually be exhausted. If L n uc 
is the typical inter-nuclear spacing, then the time to de- 
plete the ambient supersaturation, t^gp (by growth of 
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FIG. 5. Size of the stable solid droplet X2 at 
fixed time to (to ~* °°), against supersaturation a; 
a^ 1 — (<j> a — <f) g )(n / Dg) 1 ^ 2 . For a 1 < a c , the interfaces curva- 
ture-recombine at early times, and conventional growth takes 
place, at a rate proportional to a. For <r' > <r c , growth of the 
stable phase is increasingly suppressed. (Note that a c — » at 
the triple point.) 

either solid or metastable liquid) is of order L nuc /D g . 
The subsequent time-scale for metastable regions to col- 
lapse and be replaced by crystals is of order (c'/o'c) 2 t c j e p- 
This may explain, for example, the suppression of rapid 
crystallization observed beyond a sharp boundary jnj in 
colloid-polymer phase diagrams H . As a' approaches a c 
from below, the lifetime of the metastable phase diverges 
roughly as 85' 1_<T 1 , consistent with a sharp onset 
of arrested behavior. 

Notice finally that metastable phases are routinely ig- 
nored in "diffusion path analysis" (DPA), a method for 
predicting the evolution of interfacial positions in the dif- 
fusive (long-time) regime by looking solely at the equilib- 
rium phase diagram. DPA has been widely applied to 
systems undergoing dissolution pOfl . We find that, under 
dissolution conditions (ja < 0), the g/s interface is al- 
ways diffusively bound ]l4| ] . This may account in part for 
the success of DPA in that context. 

In summary we have shown that a metastable phase 
of intermediate density can radically influence the dy- 
namics of phase ordering in systems with a single con- 
served order parameter. Such an intermediate can ki- 
netically arrest the formation of a dense phase. A fuller 
treatment of these phenomena will of course require a 
three-dimensional approach (involving explicit treatment 
of contact lines and interfacial curvature). It may also 
require inclusion of additional order parameters to de- 
scribe (say) local crystallinity (implicitly treated above as 
rapidly relaxing). We leave these issues for future work. 
However, we believe that our one-dimensional treatment 
offers a strong physical picture of the arrest of domain 



growth in conserved systems caused by metastable states, 
and that this will remain an important part of any com- 
plete theory. 
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